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Abstract 

Although the weak nonleptonic amplitudes of the Standard Model are 
notoriously difficult to calculate, we have produced a modified weak 
matrix element which can be analyzed using reliable methods. This 
hypothetical nonleptonic matrix element is expressible in terms of the 
isovector vector and axialvector spectral functions pv(s) and pa(s), 
which can be determined in terms of data from tau lepton decay and 
e + e~ annihilation. Chiral symmetry and the operator product expan- 
sion are used to constrain the spectral functions respectively in the 
low energy and the high energy limits. The magnitude of the ma- 
trix element thus determined is compared with its 'vacuum saturation' 
estimate, and in the future may be accessible with lattice calculations. 



1 Introduction 



To determine the origin of the AI = 1/2 rule remains one of the unsolved 
problems of the Standard Model. Despite years of effort involving the- 
oretical techniques ranging from quark model calculations to lattice-gauge 
studies, a quantitative understanding of the AI = 1/2 enhancement is still 
lacking. The reason for this can be seen from the structure of the AS = 1 
nonleptonic weak hamiltonian, 

n*£ =1 = fv c * d V us J d 4 x iD^x, M W ) T(jQ i2 (x) JI +15 (0)) , (1) 

where iD^ u (x, Jl%) is the VT-boson propagator, g2 is the gauge coupling 
strength for weak isospin and j£ are the left-handed quark weak currents, 

^ = ^(1 + 75)9 , (2) 

where q = (u d s). Strong interaction effects are present in Eq. (|]) at all 
energy scales up to the W-boson mass. It is this fact which signals the 
difficulty of dynamically reproducing the AI =1/2 rule. 

However, experience from both experiment and theory has taught us 
much about low-energy QCD. It is possible that a hybrid approach which 
combines data with theory might provide insight as to the source of the 
AI =1/2 effect. The calculation to follow will constitute a first step in that 
direction. It involves just the vector-vector part of the operator in Eq. (|l|) 
(note that we drop the KM factors) 

n ™ =1 = f / d ' X iD ^ X > Mw) T (^2(^)^5(0)) • (3) 

and its matrix elements taken between single-particle pseudoscalar meson 
states, 

M K -^~ = (7r~\nt s=1 \K~) , 

Mgo^o = (A-Hi S=1 \K°) . (4) 

Our goal is to determine the matrix elements Mg n as accurately as possible, 
despite the complications of the strong interactions, 

The main theoretical tool will be chiral symmetry. It is easy to show that 
the single-meson matrix elements of the chiral operator of Eq. ([!]) vanish in 
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the soft meson limit. By contrast, the Mk-w of Eq. (0) have a nontrivial 
form in the chiral limit because they involve a vector-vector operator. @ By 
working in the soft meson limit, we can utilize a spectral representation to 
study the matrix elements. Constraints on the low-energy and high-energy 
limits of the spectral integrals are obtainable from chiral symmetry and 
QCD sum rule methods, and data can be used to fill in much of the rest. It 
is in this sense that the 'anatomy' of these matrix elements are revealed. 

Although the Mr-k are thus a kind of theoretical laboratory from which 
we can learn something of value, they themselves do not reproduce the 
AI = 1/2 signal. This can best be seen in terms of an effective lagrangian 
description. Under the usual classification scheme of SU(3)l x SU(3)r, 
the operator of Eq. (|3|) transforms as a component of an (8, 8) family of 
(parity-conserving) operators. HI An SU(3) invariant lagrangian which can 
characterize the momentum independent parts of the transition amplitudes 
is 

£ff = gTt (A a £A 6 £t + A 6 £A a £t) , (a, 6 = 1, . . . 8) , (5) 
where g is a constant, £ represents the meson fields, 

£ = exp {i^cK/F) , (c = 1, . . . 8) , (6) 

and F is the meson decay constant. For the nonleptonic weak (8, 8) opera- 
tor, the ratio of AI = 1/2 and AI = 3/2 amplitudes is fixed. Thus, given 
the isospin decomposition 

M K o^ w o = -^M 1/2 + V2M 3/2 , (7) 



it is straightforward to show 



M 3/2 _ 2 



2 Derivation of a New Chiral Sum Rule 

In the following, we shall express the K ^ ir transition amplitudes in the 
form of a spectral representation. Throughout, all our calculations will be 
performed in the chiral limit of massless u, d, s quarks. The first step in the 
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process is to invoke the current algebra results o 
limfr- \T (V^ l2 (x)v: +t5 (0)) \K~) 



~(0\T (^T(x)T^(O) - A£(z)AZ(0)) |0) , (9) 



and 



lim^lr^^x)^^)) \K° p ) 
3 



V2F 2 



(0|T(lf(x)^(0)-^(x)^(0))|0) , (10) 



where we take the vacuum to be an SU(3) singlet. Note the dependence 
upon the difference of vector and axialvector terms. These quantities can 
be expressed in terms of spin-one spectral functions pv,a( s ) as defined byl 

{0\T{V£(z)V?{0))\0) = 

d 4 p f-W-* 



poo p ( j'i n p—iyx 

iS ab / ds p y (s) (-a<r - 9^) / 6 (11) 

Jo J (27T 4 p z - s + ie 



and 



(0\T(A%(x)AU0)) 10) = -i6 Qb FZd»ff' 



(2vr) 4 p 2 
d 4 p e - *-* 



+^ ab / ds p A (s)(-s^ " ^5") / TTTZTI 2 „ , • (12) 



(2tt) 4 p 2 + ie 
(2vr) 4 p2-s + ie 

Note that in the chiral limit, the spin contribution is given entirely by the 
pion pole. It is thus possible to write the M. ^ n in spectral form, 

M *-~ - = iiTK^- 4 and = ~w A - (13) 



where 



A = M^I dssHn[-^) Ms \- pA{s) . (14) 



The combination of vector and axialvector spectral functions in Eq. ( |i~4| ) 
is reminiscent of similar forms appearing in other well-known sum rules viz. , 

°° ds M«)-pa(«) = _ 4Zio (15) 
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ds (pv(s) - pa(s)) 
ds s (pv(s) - pa(s)) 



ds sin f ) (py(s) - p A (s)) 



Fl , (16) 
, (17) 

V()ir2F hrnl ± -ml») ■ (18) 



3e 2 



In the first sum rule, Liq is related to the renormalized coefficient L^(p) 
of an 0{E 4l ) operator in the effective chiral lagrangian of QCD§, 



z 10 = 4oV) + 



144 



7T- 



ln 



mt 



+ 1 



-6.84 x 10 



-3 



(19) 



The next two relations are respectively the first and second Weinberg sum 
rules. M Finally, Eq. (|l^ ) is the formula for the vr^-vr mass splitting which 
was derived long ago using soft-pion methods. & Although containing an 
arbitrary energy scale A, this last expression is actually independent of A 
by virtue of the second Weinberg sum rule. As one would expect, Eq. ( |i~4| ) 



reduces to the structure of Eq. (18) upon dropping the prefactor of and 



taking the limit — > 0. Although it is the expression in Eq. (14) which 
is central to our analysis, we shall see that these other sum rules serve as 
crucial checks on the reliability of our determination. 



3 The Spectral Functions 

As mentioned earlier, the difficulty in calculating the weak nonleptonic ma- 
trix element lies in the need to understand all the physics between s = and 
s = Myj. Fortunately, with the particular combination found in Eq. ( |l~4| ) 
(i.e. pv — pa), we are able to overcome this hurdle. We have recently pro- 
vided a detailed phenomenological overview© of p\-pA arid of the chiral 
sum rules in Eqs. (15-[l8|). Briefly, the important ingredients are as follows. 
At very low energy, the difference of spectral functions is uniquely de- 



termined by chiral symmetry to be 



[10| 



1 / Am 2 \ 

PVOO - PAW ~ 4^2 ( 1 " -f J °( s - 4m D + °(P 2 ) ■ (2°) 

At somewhat higher energies, the vector spectral function py may be deter- 
mined from both e + e~ annihilations and decay of the tau lepton. As a 
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consistency check, the two data sets have_heen shown to give compatible in- 
formation in the p(770) resonance region. 113 The data reveals that as energy 



is increased above threshold, py is dominated first by 2ir and then 4tt reso- 
nances. At even higher energy, multipion production leads to a continuum 
component which ultimately approaches the asymptotic QCD prediction, 



For the axialvector spectral function pa, things proceed similarly, except 
that one must rely solely upon tau decay data. Also, above threshold it is 
the 3tt sector which first occurs and the dominant resonance contribution is 
that of ai(1260). In the large energy limit, QCD predicts that pa = Pv to all 
orders in perturbation theory. The difference between py an d pa arises 
from non-perturbative effects and may be estimated by using the operator 
product expansion. The asymptotic energy dependence behaves as s -3 , 
and in the approximation of vacuum saturation we find 

8 as(y/aSqq}l „ 3.4 x IP' 5 GeV 6 
py{s) - pa(s) ~ ^ ~ ~ 3 , (22) 

where we take a s (5 GeV 2 ) ~ 0.2. In view of the tiny numerator, this 
asymptotic tail provides only a very small contribution to the sum rule 
integrals and thus its precise value is not very important. 

A phenomenological analysis of the chiral sum rules has been performed 
in Ref. |J. Both numerical and analytical approaches have been employed 
to convert the empirical knowledge of py — pa into statements about the 
sum rules. Figure 1 displays a typical solution found there for py - pa- It 
represents a fit of the r — > 2ir, 3ir, 4-7T decay spectra and branching ratios and 
e + e~ — > 2tt, 4tt cross sections. Other solutions are found to have a highly 
similar appearance, and we refer the reader to Ref. || for details. As a 
whole, the solution set has the correct high energy and low energy behavior, 
as well as reproducing the correct values of the four chiral sum rules listed 
in Eqs. (17)-(20). Of course, the phenomenology is based on existing data 
with attendant uncertainties, especially with respect to the precise values of 
the 3ir and 4-7T branching ratios in tau decay. The form of the curve in Fig. 1 
is expected to be subject to minor refinements as future experimental results 
become available. However, such a determination of the spectral functions 
allows us to calculate the weak matrix element A, and we ascertain its value 
to be 

A = -0.062 ±0.017 GeV 6 . (23) 
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The error bars we assign about the central value is our estimate of the 
uncertainties associated with the current data base. We find that the most 
important contributions to the dispersive integral in Eq. (|l4j) arise from the 
'intermediate' energy region, 1 < s(GeV 2 ) < 10. This is seen to be a result 
of the factor of s 2 in the integrand which suppresses the effects at low s and 
of the vanishing of py - pa as s — > oo which controls the high energy region. 

We conclude this section with a model-dependent estimate of the ampli- 
tude A. Since rigorous calculation of nonleptonic matrix elements has tra- 
ditionally not been available, various approximation schemes have instead 
been employed, among them quark models and the vacuum saturation ap- 
proach. For definiteness, we shall adopt the latter approach here. The short 
distance expansion 

T(vf_ a (x)Vfa(0))=V£^(0)Vfa(0) + 0(x) , (24) 

along with the evaluation 

A^(x,M w ) = ^ , (25) 



allows us to write 



M R „~-^(K\d il »u i u j ^s j \K) . (26) 



In the vacuum saturation approximation, Eq. (26) is determined by perform- 
ing a Fierz transformation and inserting the vacuum intermediate states. 
Remembering that we are working in the chiral limit, we find 

Aac = -^<««>o - -0-033 GeV 6 . (27) 
This is about half the spectral evaluation given above. 



4 Concluding Remarks 

To our knowledge, this is the first time that anyone has undertaken evalua- 
tion of a weak nonleptonic matrix element using a battery of methods which, 
at least in principle, are both theoretically sound and practical to apply. The 
resulting analytic control allows us to be able to identify the most important 
of the underlying physics. We have found that the dominant contributions 
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to the matrix element come from the region of intermediate energies, which 
is the one over which theorists have the least control. That is, at the low- 
est energies, rigorous predictions can be made using chiral symmetry, while 
at the highest energies methods of perturbative and nonperturbative QCD 
may be invoked. However, the intermediate energy region continues to resist 
the development of reliable analytic methods. Fortunately, the existence of 
a data base allows this range to be estimated. We have done the best that 
can be accomplished at this time. As the quality and quantity of data im- 
proves, we expect to be able to reduce accordingly the uncertainty in our 
determination. 

It will be interesting to compare our evaluation of with those 

of computer simulations in lattice gauge theory. It could well be that our 
phenomenological determination, even with the presence of the aforemen- 
tioned uncertainties, is still more accurate than those of present day lattice 
studies. The energy range 1 < s(GeV 2 ) < 10 is a problematic one for lattice 
simulations because the energy scale is comparable to the inverse lattice size 
used in current computations. Lattice artifacts may then be significant. In 
addition, although today's lattice studies are becoming more accurate in the 
prediction of resonance masses and couplings, they are not yet competitive 
in matrix element analysis with the experimental determinations which were 
the prime ingredient of our evaluation. 

The weak matrix element that we have studied does not immediately 
provide the key to solving the riddle of the AI = 1/2 rule. However, it 
does provide an interesting theoretical laboratory for the development of 
improved calculational methods which may prove to be efficient in addressing 
this obdurate puzzle. Additional work on this approach is underway and will 
be reported on in future publications. 

The research described in this paper was supported in part by the Na- 
tional Science Foundation. 
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Figure Caption 



Fig. 1 The spectral function py — px 
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